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Abstract 

We present joint threshold and recoil resummed transverse momentum distribu- 
tions for heavy quark hadroproduction, at next-to-leading logarithmic accuracy. We 
study the dependence of these distributions on the production channel, the color 
configurations and the differences with the pure threshold-resummed distribution. 



1 Introduction 



The formalism PjJ |2] of hadronic cross sections for the joint resummation of distributions 
singular at partonic threshold and at zero recoil has so far been applied to only a few 
processes. The most recent studies involve processes that proceed at lowest order through 
a 2 — > 1 electroweak (Z/W production [3 ) or Yukawa interaction (Higgs production 0j). 
For these cases, the observables are the production cross sections at fixed mass Q and 
measured Qt- Partonic threshold is defined by z = Q 2 /s = 1, where s is the partonic 
center of mass energy squared, and zero recoil by Qt = 0. At any finite order, the distri- 
butions take the form of plus-distributions [ln fc (l — z)/(l — z)] and \\xi k (Q / Qt) / Qt\ 
In these observables the latter distributions enter in the physical cross sections, whereas 
the former are defined, after factorization, in the context of a perturbative analysis of the 
hard scattering. 

The basis for the joint threshold and recoil resummation for these processes and those 
proceeding at lowest order through a 2 — > 2 (QCD) reaction was given in Ref. j2J. The first 
application of the formalism to the specific case of prompt photon hadroproduction was 
presented in Ref. There the photon transverse momentum spectrum was analyzed, 
and a preliminary numerical study performed. Recoil corrections were found to have a 
significant impact on the photon spectrum, both at large and at small values of pt- For 
such 2^2 processes, the formalism implements the notion that, in the presence of QCD 
radiation, the actual transverse momentum produced by the hard collision is not but 
rather pt — Qt/^i with Qt the total transverse momentum of soft recoiling partons. In 
fact, the joint-resummed partonic pt spectrum has the form of a hard scattering cross 
section as a function of p' T = \pt — Qt/2\, convoluted with a perturbative, albeit resummed 
Qt distribution. The extreme situation Qt = 2px in which all transverse momentum is 
produced through recoil leads to a singularity in the hard scattering function, equivalent 
to the singularity at zero pt in the prompt-photon Born cross section. We note however 
that a recently proposed extension [6j of joint resummation avoids this singularity. 

In this paper we apply the joint resummation formalism as given in Ref. [2] to another 
prominent 2 — > 2 scattering observable: the pt distribution of heavy quarks produced in 
hadronic collisions. Key differences with the prompt-photon case are, first, the presence 
of the heavy quark mass m, preventing a singularity in the hard scattering function when 
Qt = %Pt and, second, the possibility of multiple colored states for the produced top 
quark pair. We derive joint-resummed and threshold-resummed transverse momentum 
distributions at next-to-leading logarithmic (NLL) accuracy, and study for different pro- 
duction channels the dependence on color configurations, as well as differences among the 
two resummations and the exact next-to-leading order (NLO) distribution. 

Other recent studies for heavy quark pt distributions based on finite order expansions 
of threshold resummation can be found in Refs. [3 |S] . 

Relevant formulas are derived in section 2. In section 3 we present numerical studies 
of the dependence of the joint resummation effects on the production channel, on color 
configurations and on the heavy quark flavor. We also assess the differences with the pure 
threshold-resummed distribution. Section 4 contains our conclusions. 
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2 The observable 



We consider the inclusive pt distribution of a heavy quark produced via the strong inter- 
action in a hadron-hadron collision at center of mass (cm) energy 

h A {p A ) + h B (p B ) -> Q(p c ) + X , (1) 

where h A ,B refers to the two incoming hadrons, Q to the detected heavy quark and X to 
the unobserved part of the final state which includes also the heavy anti-quark Q. The 
lowest order QCD processes producing a heavy quark are 

s(po) + q(pb) -> Q(Pc) + Q(Pd) , ^ 
g(pa) + g(Pb) -> Q(p c ) + Q(Pd) , 

at cm energy y/§ = y/^bS with £ a ,&> parton momentum fractions. Exact higher order 
corrections to the differential cross sections for these partonic processes have been com- 
puted to NLO Ell HU E21- To any order (T3] the observable may be written in the 
following factorized form (up to power corrections behaving as l/p T ) 



dp 7 



yZ I d ^ d £b <fra/ A (L, ^)<Pb/B(£b, y) d(Tab -* Q+x fc g> ^ a 3 (n),pr) , (3) 
~~r Jo a Pr 



with da a b^Q+x/dpT the partonic differential cross-section, 4> a / A and <fib/B parton densi- 
ties, and /i the factorization and renormalization scale. The purpose of this study is to 
investigate certain effects of soft gluons in these higher order corrections for the heavy 
quark p T distribution. 

The first of these are threshold enhancements, which essentially involve the energy 
of soft gluons. In the context of the factorization @ we define hadronic and partonic 
threshold by the conditions S = 4m^ and s = 4m^, respectively, with rriT the transverse 
mass a/ 'm 2 + p\. It is convenient to define the scaling variables 

x t = 7T~ ' x t = c c c ' v*/ 

so that hadronic (partonic) threshold is at = 1). The higher order corrections 

to the partonic cross section da a bjdpT contain distributions that are singular at partonic 
threshold. Threshold resummation organizes such distributions to all orders. 

The second of these are recoil effects, resulting from radiation of soft gluons from 
initial-state partons. We wish to treat these effects in the context of joint threshold and 
recoil resummation. To this end we employ the refactorization analysis of Refs. E], 
which enables us to identify a hard scattering with reduced cm energy squared Q 2 and at 
transverse momentum Qt with respect to the hadronic cm system. This hard scattering 
produces a heavy quark with transverse momentum 

Pt = Pt~^- (5) 

The kinematically allowed range for the invariant mass Q of the heavy quark pair in this 
hard scattering is limited from below by 2m' T = 2 a/ m 2 + p'£ so that threshold in the 
context of joint resummation is defined by 



A refactorization analysis on the lines of the one performed in Ref. [2j leads to the following 
expression for the observable in Eq. (J5J) 

d(JAB -^ x = f d 2 Q T dQ 2 9{JL - \Q T \) toAB-^+x (?) 
dp T J dp T dQ 2 d 2 Q T 

where 

li '~ ( - +X - ^ [ d£ad£b<Pa/A(£a,li)<f>b/B(£b,fi) 



dpTdQ 2 d 2 Qx 



ab=qq,gg 

x J dx a d 2 k a c a / a (^-,k^ j dx b d 2 k b c h/b (^ 1 k b 

x iw dWs w> + ^(^ ) 

x 1*5(1 - Q 2 /S - (1 - x a ) - (1 - x b ) - w s )6(Qt + K + h) 



is known as the profile function. Note that we restrict ourselves to relatively soft, per- 
turbative recoil contributions by integrating over Qt in Eq. (JZj) such that \Qt\ < fi- 
Implementing a cut-off p, is of course a rather crude approximation to what should be 
a consistent matching to finite order results. The implementation of such a matching 
procedure is however beyond the scope of this work. 1 We shall discuss the limitations on 
the integration variables such as Qt further in the next section. 

Let us comment on the right hand side of Eq. (jHJ, moving from bottom to top. The last 
line implements transverse momentum conservation, as well as the decomposition of the 
(normalized) above-threshold energy 1 — Q 2 /S into contributions from collinear radiation 
(1 — Xa )b ) and from wide-angle soft radiation (w s ). The definition of the observable and the 
refactorization ensure [2] that only the initial state contributes to the recoil of the hard 
scattering. The line above the last displays the hard scattering functions C$ da^ /dp' T , 
one per color state (singlet or octet) in which the heavy quark pair can be produced. In 
a diagrammatic representation all its lines are off-shell by at least an amount p' T . The 
factor dp' T /dpT accounts for the phase space difference between the hard scattering in the 
factorized expression Eq. (j3J) (da^/dpr) and here, where da^ jdp' T is appropriate. The 
functions S[ and S' s summarize the effects of coherent wide-angle soft radiation, and are 
sensitive to the color structure of the hard scattering. In particular, the effect of soft 
radiation for our observable is such that, to NLL accuracy, there is no mixing between 
singlet and octet color states, as was already observed for the case of the total heavy-quark 
cross section near threshold [14J. 

Moving up further in Eq. (jHJ), the functions c a / a are defined in Ref. and are related 
to the density of parton a in a parton of the same flavor at fixed energy fraction x a and 
transverse momentum k a . The prime on the S function indicates that factors \/Udy, 
the square root of the Drell-Yan soft function, have been introduced in equal numbers in 
numerator and denominator and then conveniently redistributed between the heavy quark 
soft functions and the c functions in order to remove gauge dependence in the individual 
functions. Finally the top line displays the convolution with standard parton distribution 



1 Note that such a matching in the extension of Ref. [H] is straightforward. 
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functions, as well as the sum over initial state parton flavors, in which the gq, gq channels 
have been neglected because they are suppressed near threshold. 

A convenient way to implement the energy and transverse momentum conservation is 
by passing to transform space, where the impact vector b is Fourier conjugate to Qt, and 
the variable N is Mellin conjugate to the threshold variables in Eqns. (j3J) and (JBJ). This 
converts the convolution in Eq. (JSJ) to a product. We aim at NLL accuracy in transform 
space, that is we account for all contributions exp{a™L n+1 } and exp{a™L n }, where L can 
be either In N or In b. 

Using their renormalization group scaling properties, moments of the soft functions 
may, to next-to-leading logarithmic (NLL) accuracy, be written as jT3 EH] 



/8 



Q_ 

Ay 



a s (ii) ) = S' 1/a (l,a B (Q/N))exp 



Q/N du' 

-^2Rer 1/8 (a s ( /U / )) 



(9) 



The solution needs careful definition in order to avoid values of N for which a s (Q/N) 
is singular. The iV-integration contour we choose indeed assures this, and anyhow, at 
NLL accuracy, we can approximate the functions Sy s (l,a s ) by their lowest order ex- 
pressions. The form of the one-loop soft anomalous dimensions Ti and Tg we need is 
straightforwardly derived from the expressions in Ref . fSl GUI El • Their real parts are 



2ReF 1 (a s ) = 2Cp (ReLg + 1) , 

71 

2Rer 8 (a s ) = 2 Re T x (a a ) + — C A In + ReL^ 

7r V m 2 



(10) 
(11) 



ReL, 



l + {3 2 



In 



m 2 /m 



(12) 



2f3 V 

We postpone a discussion of these expressions to further below. 

To NLL accuracy we may approximate both dp' T /dpx and the C$ functions by 1, so 
that integrating over Q 2 , and, with manipulations similar to those of Ref. [21 we arrive at 



dpxd 2 Q 



ab=qq,gg 



g -2 C F t(N) (Re^+1) CAt{N) ( Xn ^ +Ld )\ S 



M((N) + M£(N)e 



^S 2 V / V 4 ( m ' + \Pt-Qt/2\' 



N+l 



(13) 



where the inverse Fourier transform and Mellin transform are explicit. Notice in particular 
the last factor, which provides a kinematic link between recoil and threshold effects. The 
exponential functions E a b E] to next-to-leading logarithmic (NLL) accuracy are 

E ab (N,b)= [ Q ^L[A a (a s (fx')) + A 6 ( as (//))]21n^ - gb 2 , N = Ne» , (14) 

where the coefficients A a and are taken from Ref. We also added to the perturbative 
exponent the non-perturbative (NP) Gaussian smearing term — gb 2 . 
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To simplify notations, we have combined in Eq. (|13|) S'^Sg with the Born matrix 
elements that contribute to da^/dpx to build up the two Born functions M^(JV) and 
M|(iV), and we have introduced the evolution variable 

The Born functions Ml(N),Mg(N) are the Mellin moments of the lowest order heavy 
quark production matrix elements for either the qq or gg channel, the index labeling the 
color-state of the heavy quark pair: 

r l (r 2 ) N 



s N c T(N + 5/2) V V ; 



m 2 " 
m 2 T/ 



(16) 



Jo v 1 — xf. 



, ^ , 2T R ^r(iV + 1) / N + l 2 (rn^_rn* 
^ s N C {N 2 -1) T{N + 3/2) V N \m%, m T 



1 f£.2\N 



M* gi8 (N) ee / ^^^ 8 (f 
Jo a/1 — 



167r 2 ag y^Fr(iV+ 1) 

TV 2 — 1 r(iV + 3/2) |_V"~' ^7 V N 

C A fN 



2N + 3 V N 



2(iV+l) 



3 


m 4 \ " 






to 


4 

m 







with Mf and M| the square matrix element for the production of a hard QQ system in 
a singlet and octet color state respectively. We point out that, at this order, in the qq 
channel the heavy quark-antiquark pair is produced only in an octet state. 

In a refactorization analysis, the expressions 2Re r x and 2Re T 8 are related to the (real 
parts of) anomalous dimensions of certain operators, composed of Wilson lines jTHl EH EH] , 
as the structure of Eqs. (JUJ) reflects. Because for our case imaginary parts of the virtual 
corrections (Coulomb phases) to the anomalous dimension cancel in the cross section, the 
soft function can be given a natural probabilistic interpretation. 

From this viewpoint, at NLL accuracy, soft wide-angle gluons can be considered as 
being emitted independently from the external partons, and approaching threshold is 
equivalent to suppressing all soft gluon emissions with energies above Q/N. Furthermore, 
due to the fact that there is no mixing between different color structures, soft gluon emis- 
sion exponentiates straightforwardly. The exponential in Eq. © can then be interpreted 
as the probability of not emitting any wide-angle soft gluon with energy above Q/N. We 
now offer some comments on the explicit expression of the anomalous dimensions in some 
limiting cases. 
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• In the case in which the heavy quarks are at rest in their cm frame, i.e. (3 — > 0, 
ReLg — > — 1, and 2Rer! vanishes, corresponding to the fact that the combined heavy 
quark-anti-quark pair forms a color- singlet. On the other hand, 2Re T 8 — > — oi s CaI 7r, 
which reflects the fact that the quark- ant iquark pair produced in an octet state 
radiates as a single object rather than two separate particles. 

• For positive moderate values of (3 (0 < (3 < 1) we observe that r§ contains the term 

In ( m T \ = \ n ( (2p Q • p c ) (2p b • p c ) 
\m 2 J \ (2p a ■ p b ) m 2 

which is the standard contribution to the anomalous dimension from the qqg an- 
tenna. 

• In the extreme region 3> m ((3 — > 1), gluons can become effectively collinear to the 
heavy quarks. In this case Lg approximates lnp T /m 2 , a large collinear logarithm. 
The resulting logarithmic contributions are identical in Ti and T&, reflecting the 
fact that collinear radiation depends only on the color charge of each hard emitting 
parton, and is not sensitive to the color or geometrical structure of the event. In 
this large px limit, collinear logarithms should be resummed to all orders, and a 
different analysis is required [201 EE] • 

This completes our discussion of the joint-resummed transverse momentum distribu- 
tion. The threshold-resummed result can now easily be derived, by substituting Eq. (|13|) 
into (0) (after the performing the Q 2 integral) and neglecting Qt in the last factor in 
Eq. (|13|) . Then the Qt integral sets b to zero everywhere, yielding the threshold-resummed 
result. 




3 Top and bottom quark spectra 

In this section we demonstrate numerically the soft gluon effects for top and bottom quark 
transverse momentum distributions for collider and fixed target kinematics respectively. 
Rather than performing a detailed comparison with data, our aim is to examine the 
numerical behavior of our formulas in realistic settings. We leave therefore also to a 
future study the application of our results for the description of charm quark/charmed 
meson pt spectra in fixed target kinematics. 

Let us first comment on the allowed integration range of the recoil momentum Qt- 
In Eq. (J7J) we have simply imposed an upper limit fx on the modulus of this transverse 
momentum vector. There is however another constraint on Qt, because the recoil factor 
at the end of Eq. (JT3*j) implies the condition 



4(m 2 + \p T -Q T /2\ 2 ) 

which expresses the fact that the available hadronic energy be more than the minimum 
mass of the final state in the recoiling hard scattering. Clearly, for the maximum value 
p T = S/4 — m 2 , Q T has no phase space left. For a given value p T and (large) jl the allowed 
range can be chosen to be the overlap of the two disks Di 2 in Qt space 

D 1 : \Q T - 2p T \ 2 < S - Am 2 , D 2 : \Q T \ < Jl . (21) 
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or one can choose jl such that disk D 2 fits in Di. We chose the former option. Notice 
that due to the heavy quark mass, the expression in Eq. (JT3J) stays finite at Qt = 2pr , 
in contrast to the prompt photon production case analyzed in Refs. |2J E] • 

For top quark production we chose the cut-off /2 = 200 GeV. This might seem quite 
large as a cut-off on recoil momentum, given that the joint-resummation formalism as- 
sumes that recoil radiation is soft. This large value is however relatively innocuous, 
because, as we will see, for the joint-resummed cross section the dominant Qt values are 
quite moderate, typically a few percent of the hard scale. The benefit of this large cut-off 
value is that, while recoil effects are still dominated by soft gluons, it allows us to recover 
numerically the threshold- resummed result, in the manner described at the end of section 
121 For bottom quark production we choose p, = 30 GeV, the motivation for this large 
value being similar to that for the top quark. 

Before we show results, let us state our default choices for various input parameters. 
For these and other plots we use the NLO GRV parton density set of Ref. [22] corre- 
sponding to a s (Mz) = 0.114, with the evolution code of Ref. |22|, changing flavor number 
at = m c (1.4GeV), mb(4.5GeV). We chose the factorization and renormalization scale 
equal to the transverse mass uit- For the non-perturbative Gaussian smearing parameter 
g in Eq. (J 14)) we took g = lGeV 2 . Furthermore, to avoid spurious singular behavior at 
the subleading level, we used, following [3| 

N - bOe m 

\{bQ,N) = b+ Y+7fi/N' < 22 > 

with 7] = 1/4. 

We begin with top quark (m = 178 GeV) production, and compare in Figure ^ joint- 
resummed, threshold-resummed and the exact LO and NLO calculations for the top quark 
Pt spectrum produced in the qq channel at the Tevatron (y/S = 1.96 TeV). We observe 




Figure 1: Top quark spectra and .fT-factors for the qq channel. 

that, while the resummed and NLO curves are close for small and moderate px (the 
inset provides a somewhat better view of the low pt region), for large values the re- 
summed curves depart significantly from the NLO curve. Of course, cross sections for 
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top quark production at such large px at the Tevatron are far too small to be measured, 
so that our plots at large px have only theoretical interest. For such large px values, 
the hadronic threshold, defined in Eq. (@J), approaches the partonic one, where larger N 
values dominate, a prerequisite for seeing significant effects for both resummations. The 
enhancements relative to the Born cross section are shown in the form of a K-factor 2 in 
Figure H As expected, threshold resummation produces an overall enhancement of the 
cross section that increases with increasing px, yielding e.g. a 35% enhancement over 
NLO at px = 800 GeV. Joint resummation almost doubles that effect, which can be un- 
derstood from the arguments given in the introduction: the joint-resummed enhancement 
at large px effectively constitutes a smearing of the threshold-resummed px spectrum by 
a resummed recoil function. 



p T = 500 GeV p T = 900 GeV 




Figure 2: Top quark profile functions (Eq. (|T3|)) for the qq channel, at px = 500 and 900 
GeV. 

Before showing results for the gg channel, we exhibit in Fig. |2] the Qx profile of Eq. (|13|) 
for two rather large px values, in analogy to Fig. 1 in Ref. [Hj- In agreement with Fig. 
we observe a small enhancement over the threshold-resummed result, in particular at very 
large px- The peak of the distribution is at relatively soft values of Qx, a few percent 
of the hard scale rrix- For px = 500 GeV Qx reaches /i before the bound in Eq. (j^UJ), 
so that there is no divergence for larger Qx- For px = 900 GeV there is a region of Qx 
phase space where the bound in Eq. (f2T)J) is saturated and the cross section diverges. The 
onset of that singularity is visible in Fig. but because of our choice fi = 200 GeV it 
does contribute a small amount to the i^-factor in Fig. ^ at the largest px value. We also 
observe that, in case of top production, for px < 863.7 GeV, the disk D 2 in eq. (|2*T|) is 
fully contained in the disk Di which represents the bound in Eq. ()20|). Therefore below 
that value, which includes all px values accessible experimentally, all the enhancement 
of the joint resummed cross section with respect with the threshold resummed is to be 
ascribed to recoil effects. 

In Figs. |3]we show results analogous to Fig.Q now for the gg channel summed over the 
two color states, while in Fig. H]we show the K-factors per color state. The K-factors show 

2 Normalized to the Born cross section. 
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Figure 3: Top quark pt spectra and corresponding fT-factors for the sum of color-sing 
and color-octet gg channels. 




Figure 4: ^-factors for the incoming color-singlet and color-octet gg channels. 
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a significant enhancement for joint resummation with respect to threshold resummation, 
increasing with increasing px, and strongest for octet production. The resummed spectra 
are harder than the fixed order ones. The increased enhancement of the gg channel over 
the qq channel is due to the larger value for A g compared to A q in Eq. (fTlj) . 

We have verified that both resummed cross sections have a reduced scale dependence 
with respect to LO, as expected on general grounds flQ I25j. The reduction in the re- 
summed cross section is somewhat less than in the exact NLO cross section, suggesting 
that a matched NLL-NLO calculation would do better still in this regard. Doing such a 
calculation is however beyond the scope of our paper. 

In order to assess the relevance of the last term on the right hand side in Eq. (fTTj) . 
which encodes the differences between contributions from octet and singlet soft anomalous 
dimensions, we computed the px spectrum for the qq channel with and without this term. 
We found an effect ranging from 8% enhancement with this term included at lower px to 
about 10% at higher px- For the octet gg channel we found the enhancement over the 
singlet-only to be 12 to 15%. 

Turning now to the bottom quark px spectrum, we took as kinematic conditions, for 
illustrative purposes, pp collisions at the HERA-B cm energy (41.6 GeV). In these settings 
the gluon channel is more dominant, and effects of joint resummation are more noticeable 
than for top production in the above. Fig. El is analogous to Fig. [T] but the effect of joint 
resummation are stronger now, which holds as well for the gg channel in Fig. |H1 In Fig. 
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Figure 5: Bottom quark px spectra and K-f actors for the qq channel. 



for the qq channel a significant effect at lower p? of virtual contributions in the NLO cross 
sections, not present in the resummed versions is visible, which is much less in the gg 
channel. The .fT-factors are large (which is not unusual El Ej for this kinematical 
configuration), with a clear distinction between threshold and joint resummation. The 
Qt profile of the enhancement at two values of pt is shown in Fig. [7| Here is there no 
divergence, because the bound Eq. (|2*D*|) is simply not reached. 
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Figure 6: Bottom quark spectra and i^-factors for the gg channel. 




Figure 7: Bottom quark profile functions (Eq. (fT3|) ) for the gg channel, at pt = 5 and 8 
GeV. 
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4 Conclusions 



In this paper we derived the joint threshold and recoil resummed heavy quark transverse 
momentum distributions, to NLL accuracy. As a corollary we obtained the pure threshold 
resummed result as well. The resummed distributions contain two mutually incoherent 
components, associated with different color states (singlet and octet) of the heavy quark 
pair. These components differ by their underlying Born process and by anomalous di- 
mensions governing the pattern of wide-angle soft emission. We provided a probabilistic 
interpretation of these anomalous dimensions. 

For top quark production at the Tevatron, and bottom quark production at HERA-B, 
we studied the dependence of these results on the color state, and the difference with 
threshold-resummation. We found that joint resummation gives results noticeably differ- 
ent from threshold resummation, and that for a realistic range of transverse momentum 
the kinematic singularity present in the formalism, which is very visible in applications 
to massless (e.g. prompt photon) particle px spectra, is well-screened by the heavy quark 
mass. Not surprisingly, resummation effects are larger in the gg channel than the qq chan- 
nel, and more noticeable for bottom-quark spectra at HERA-B than top-quark spectra at 
the Tevatron. The effect of the pure octet part of the anomalous dimension is noticeable. 
We hope that this study, taken together with recent other ones IH E] contributes to 
developing joint resummation as a viable framework for making QCD predictions. 
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